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Spontaneous  formation  of  composition-modulated  structures  is  a common  phenomenon 
in  III-V  and  II- VI  semiconductor  alloys  [1].  There  exist  two  distinct  possibilities  for  the 
formation  of  such  a structure.  Firstly,  equilihrium  domain  structures  can  be  formed 
in  closed  systems.  This  is  realized  by  long-time  growth  interruption  or  by  post-growth 
annealing.  Thermodynamics  can  be  applied  to  describe  equilibrium  structures  which 
meet  the  condition  of  the  Helmholtz  free  energy  minimum.  Secondly,  non-equilibrium 
structures  can  be  formed  in  open  systems.  These  structures  are  governed  by  growth 
kinetics. 

Although  most  of  observed  composition-modulated  structures  in  semiconductor  al- 
loys correspond  to  as-grown  samples,  there  is  a possibility  for  the  formation  of  modulated 
structures  via  annealing.  E.  g.,  experiments  on  non-stoichiometric  As-rich  GaAs  show 
the  enhancement  of  bulk  diffusion  of  As  atoms  under  annealing  and  the  growth  of 
nanometer-size  As  clusters  [2],  These  data  make  the  problem  of  equilibrium  modulated 
structures  in  semiconductors,  which  can  be  obtained  under  annealing,  to  be  a problem 
of  both  experimental  and  theoretical  interest. 

Here  we  study  the  thermodynamic  stability  of  an  alloy  Aj  _fB,  C in  the  epitaxial  film 
against  fluctuations  of  alloy  composition.  We  consider  a film  on  the  (OOl)-substrate  of 
a cubic  crystal  and  we  focus  on  the  situation  where  the  homogeneous  alloy  with  the 
composition  c = c is  lattice-matched  to  the  substrate.  Both  the  substrate  and  the 
epitaxial  film  are  elastically-anisotropic  cubic  crystals  having  equal  elastic  modulus 
tensors. 

The  free  energy  of  an  inhomogeneous  alloy  is  a sum  of  the  chemical  and  elastic 
contributions,  Finh  = Fchem+Ee h where  Achem  = / d3rfihem(c(x,y,z)),  and 


Fe  1 = 


dkxdky 


(27 


dz  / dz'  A c (kx , ky\  z)Be\  (kx  ,ky\z,z)  A c(kx , ky\  zf ) . 


(i) 


Here  A c(kx,  ky\z)  is  the  Fourier  transform  of  the  composition  fluctuation  A c(x,y,  z)  = 
c(x,y,  z)  — c,  z = 0 is  the  interface  plane,  z = h is  the  planar  stress-free  surface,  elastic 
properties  of  the  crystal  are  described  through  the  kernel  Be\(kx,  ky\  z,  z')  defined  in  [3], 
and  the  gradient  energy  (~  (Vc)2)  [4,5]  is  omitted,  since  it  is  not  important  for  alloy 
decomposition  in  the  film  [3].  The  problem  of  finding  the  equilibrium  profile  of  alloy 
composition  consists  of  two  parts.  The  first  part  is  to  find  the  criterion  that  the  ho- 
mogeneous alloy  is  absolutely  unstable,  i.e.  unstable  against  infinitesimal  fluctuations 
of  composition.  The  second  part  is  to  obtain  the  final  equilibrium  structure  of  the 
decomposing  alloy.  The  first  problem  was  solved  in  [3],  and  the  criterion  of  absolute 
instability  is 


(a2/chem  (c;r)  ) +Amin<() 


(2) 
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Here  A™n  =min  Xo(kx,  ky),  and  Ao (kx,kY)  is  the  minimum  eigenvalue  of  the  operator 

kx,ky 

Be i for  a given  k = ( kx,ky ).  The  minimum  eigenvalue  \<)(kx,  ky)  is  governed  by  the 
anisotropy  of  the  elastic  modulus  tensor  and  is  attained  for  wave  vectors  along  the 
elastically  soft  direction  [100]  or  [010].  The  value  A™nis  attained  asymptotically  as 
kh  -7-  oo.  Fluctuations  of  composition  corresponding  to  the  “soft  mode”  are  localized 
at  the  surface  z = h and  decay  away  from  the  surface, 

AcSoft  modeC^j  y<>  z|k)  ~ exp[-\kx\(h  —z)\  exp (ikxx)  , \kx\h  > 1 , ky  = 0,  (3) 


and  similar  expression  is  valid  for  the  “soft  mode”  with  k = (0,  ky).  At  high  tempera- 
tures, the  left  hand  side  of  Eq.  (2)  is  positive  for  all  compositions.  The  temperature  Tc 
and  the  composition  cq  where  the  l.h.s.  of  Eq.  (2)  vanishes  for  the  first  time  correspond 
to  the  critical  point  of  alloy  instability  in  the  epitaxial  film. 

The  final  structure  of  decomposing  alloy  was  found  numerically  in  [6],  where  the 
model  regular  solution  approximation  [7]  was  used  for  the  chemical  free  energy  of  an 
alloy.  An  equilibrium  structure  was  found,  but  neither  the  miscibility  curve,  nor  the 
dependence  of  the  period  on  temperature  were  obtained  in  [6]. 

Here,  to  obtain  th e final  equilibrium  structure  of  the  decomposing  alloy  we  focus 
on  the  vicinity  of  the  critical  point,  i.e.,  on  T ss  Tc  and  c « co,  and  use  the  Ginzburg- 
Landau  type  of  expansion  of  Fmu.  We  take  the  composition  fluctuation  is  a linear 
combination  of  soft  modes  of  the  form  of  Eq.  (3).  Since  numerical  calculations  show 
composition  profile  to  be  one-dimensional,  we  focus  on  composition  modulation  in  the 
x direction.  By  taking  into  account  the  asymptotic  behavior  of  the  eigenvalue  of  Bt\ 
for  soft  modes  [3],  Xo(kx,0)  = A™11  + B\  exp(— 2\kx\h),  one  obtains  the  expansion  for 
A F = Tjnh  - Thom  as  follows, 


A F 


’ j^dzjdx { J (c(r)  - c0)2  - ^r(c(r)  -c0)2  + ^D(c( r) 

1 ( <92/chem  ' 


2 V dc2 


(c-c0)2  - ^ar(c-co)2  + ^D(c  - c0)4 


+ 1 f_dkx 


V (2tt)./o 


dz  [A™11  + B\  exp  ( — 2 1 A: v | A ) ] |A c(kx,  0;z)|2 


(4) 


Here  r = (Tc  - T)/Tc  < 1 , a = -Tc(d3fcbem/dTd2c) , D = l(d\fchem/dc4).  We 
seek  the  composition  profile  as  a periodic  function  A c(x,y,z)  = \j2a./ D^frpfx,  z ). 
Here  <p  is  combination  of  soft  modes  of  Eq.  (3)  corresponding  to  composition  modu- 
lation in  x-direction  with  the  period  dx.  By  substituting  A c(x,y,z)  of  the  above  form 
into  the  free  energy  of  Eq.  (4),  one  obtains 


A F = Ah 


6a2r2 

D 


1 B\  exp(-2 kf]h) 

kx^h  aT(2kx>>  h) 


(5) 


where  A is  the  substrate  area,  and 
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Here  77  is  proportional  to  the  deviation  of  the  average  alloy  composition  c from  the 
critical  composition  cq,  77  = (c  — co)/v/r;  a\  is  amplitude  of  the  first  Fourier  harmonic 
of  ip.  It  follows  from  Eq.  (6)  that  J does  not  depend  on  ki°\  Minimization  of  A F from 
Eq.  (5)  with  respect  to  k{0>  shows  that  the  second  term  in  Eq.  (5)  is  small  compared  to 
the  first  one.  It  allows  to  split  the  problem  into  two  parts,  and  to  obtain  the  composition 
profile  within  one  period  of  the  structure  first,  and  to  find  the  period  afterwards.  This 
is  the  key  advantage  of  using  Landau-Ginzburg  expansion  of  Fin),  which  enables  us 
to  proceed  as  compared  to  [6].  Similar  advantage  of  using  Landau-Ginzburg  theory 
was  emphasized  in  [8]  for  a related  problem  of  equilibrium  composition-modulated 
structures  in  bulk  samples  of  quaternary  alloys. 

The  lowering  of  the  total  free  energy  of  the  system  due  to  fluctuations  of  composition, 
A F < 0,  means  that  the  minimization  of  J from  Eq.  (6)  yields  a non-trivial  solution 
providing  J < 0.  If  | ry  | < 1 / x/3,  the  quadratic  coefficient  in  Eq.  (4)  is  negative,  i.e.  the 
alloy  is  absolutely  unstable.  Then  there  obviously  exists  a solution  which  gives  J < 0. 
The  equation  77  = ± 1 / \/3  yields  the  spinodal  curve  of  the  alloy  epitaxial  film  near  the 
critical  temperature  Tc.  Due  to  cubic  terms  in  J,  a solution  providing  J < 0 may  exist 
also  in  a certain  region  of  77  where  |?7|  > 1 / x/3 - The  value  (77!  =77*,  where  the  solution 
providing  J < 0 ceases  to  exist,  corresponds  to  the  miscibility  curve.  If  I77I  > 77*, 
the  homogeneous  alloy  epitaxial  film  is  absolutely  stable.  If  l/v^3  < 77  < 77*,  the 
homogeneous  alloy  epitaxial  film  is  metastable.  Numerical  minimization  of  J yields  the 
final  structure  of  the  alloy  within  one  period.  The  form  of  the  structure  is  close  to  one 
obtained  in  shcha.  The  modulation  amplitude  is  maximum  at  the  surface  and  decays 
away  from  the  surface.  Analysis  of  J versus  77  yields  77*  = 0.784,  i.e.  it  gives  the 
miscibility  curve  in  the  vicinity  of  the  critical  point, 

7miscibi1ity(d)  = Tc  - - C0)2  • (7) 


Then,  given  J,  the  criterion  of  the  minimum  with  respect  to  kl(>)h  reads 

exp(2^0)/;)  _ B\  2 
(A'(°)/7)  = ar(— J)  aX  ' 


(8) 


To  calculate  the  period  d,  we  take,  as  an  example,  GaAsi_cSbc  where  Tc  = 740  K [3]. 
Evaluation  of  d from  Eq.  (8)  yields,  for  Tc  T = 10  K,  d = 1.4/?,  for  Tc  — T=  100  K, 
d = 2. Oh.  Thus,  the  period  is  a weak  logarithmic  function  of  (Tc  — T)  and  is  of  the 
order  of  the  film  thickness  h. 

With  further  lowering  of  T,  all  modes  of  fluctuations  contribute  to  the  equilibrium 
structure.  To  estimate  the  miscibility  curve  at  arbitrary  T,  we  use  the  ansatz  of  [9], 
where  composition  is  independent  of  z,  Ac(x,y,  z)  = c\(x).  If  we  construct  a curve 
which  bounds  the  region  on  T — c phase  diagram  corresponding  to  alloys  unstable 
against  such  a fluctuation,  this  curve  will  lie  inside  the  exact  miscibility  curve.  The 
equilibrium  phase  diagram  of  alloy  epitaxial  film  is  shown  in  Fig.  1. 

To  conclude,  we  have  solved  a non-linear  problem  for  the  final  equilibrium  structure 
of  phase-separating  alloy  in  an  epitaxial  film  in  the  vicinity  of  the  critical  point.  The 
exact  phase  diagram  in  variables  “temperature  - average  composition”  is  constructed, 
containing  the  region  of  absolutely  unstable,  metastable,  and  stable  alloys.  The  period  of 
equilibrium  structure  is  found.  An  equilibrium  composition-modulated  structure  with  a 
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Fig  1.  Equilibrium  phase  diagram  of  an  alloy  in  the  epitaxial  film.  (c0;  Tc)  is  the  critical  point  of 
alloy  instability  against  spinodal  decomposition.  1 is  the  spinodal  curve;  the  solid  part  of  the  curve 
2 is  the  calculated  miscibility  curve  (Eq.  (7));  the  dashed  part  of  the  curve  2 is  the  schematic 
miscibility  curve  for  arbitrary  temperatures;  3 is  the  miscibility  curve  calculated  by  the  ansatz  of 
Glas  [9],  (I)  is  the  region  of  the  absolutely  unstable  alloys,  (II)  is  the  region  of  metastable  alloys, 
and  (III)  corresponds  to  stable  alloys. 


one-dimensional  periodicity  in  the  surface  plane  is  a spontaneously  formed  quantum- 
wire  superlattice. 
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